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Abstract 

The concept of topological persistence, introduced recently in computational topology, finds appli- 
cations in studying a map in relation to the topology of its domain. Since its introduction, it has been 
extended and generalized in various directions. However, no attempt has been made so far to extend 
the concept of topological persistence to a generalization of 'maps' such as cocycles which are discrete 
counterparts of closed differential degree one-forms, a well known concept in differential geometry. We 
define a notion of topological persistence for 1-cocycles in this paper and show how to compute its rele- 
vant numbers. It turns out that, instead of the standard persistence, one of its variants which we call level 
persistence can be leveraged for this purpose. It is worth mentioning that 1-cocycles appear in practice 
such as in data ranking or in discrete vector fields. 
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1 Introduction 



The persistent homology, introduced recendy |f8| in computational topology, studies the longevity of topo- 
logical features as one sweeps a topological space with growing sublevel sets of a real-valued function. 
Since its introduction, the concept has been generalized in many directions, including its study for homology 
groups under different coefficient rings lfT4l . its stability under function perturbations lH, its multidimen- 
sional version 131, and other variants llIllSHTlfTOl. 

In a discrete setting, which is prevalent in practice, the topological space in question is usually taken 
as the one defined by a simplicial complex X whereas the function in question is taken as a linear map 
/ : X — )• M whose restriction to each simplex is linear. Such a pair (X, /) provides what is called a 0-cochain 
in algebraic topology. A 0-cochain can be specified by providing only a n-dimensional vector of /-values 
on n vertices of X. Such data appear frequently in applications involving shape and data analysis. Similarly, 
a more general input may specify real values on edges of the simplicial complex by presenting a skew 
symmetric n x n matrix for a simplicial complex X with n vertices. Such a data, in reverence to the term 
0-cochain, is called a 1-cochain. An interesting situation, which occurs in practice, is when the 1-cochain 
is a 1-cocycle (definition in section 2). For example, in data ranking, pairwise orderings which are locally 
consistent represent a 1-cocycle |[T5l[T6l . A sampling of a curl-free vector field on a Riemannian manifold 
provides interesting examples of 1-cocycles. 

It is natural to study the qualitative complexity of 0-cochains and 1-cocycles in analogy with Morse 
theory and Novikov theory |[T2l respectively. In particular, one would like to relate the dynamics of the 
0-cochain and 1-cocycle with the topology of X. For 0-cochains, this is precisely achieved by studying 
persistent topology. Although the concept has been generalized in many directions as we mentioned, it has 
not yet been extended to 1-cocycles. It is appropriate to mention that persistence and circle valued maps 
are also in the attention of ifTSl where the existence of a circle valued map is established using standard 
persistence. 

In this paper, we provide a definition of persistence for almost integral 1-cocycle^ defined on a simpli- 
cial cornplex and an algorithm to compute relevant numbers. First, we show that a 1-cocycle when almost 
integral □ can be interpreted as a S^-valued (circle valued) map. 

The standard persistence E] [HI can not be applied to circle valued map since sublevel sets are not 
defined. We observe that level persistence originally introduced as interval persistence in |7J can be lever- 
aged to define a notion of persistence for circle valued continuous maps. This persistence uses a lifting of 
the circle valued map which becomes a M— valued map on a lifted space of X which, in topologists' ter- 
minology, is the infinite cyclic cover associated with the cohomology class of the almost integral cocycle. 
We show how one can calculate the relevant numbers for level persistence for this lifted map by adapting 
the algorithm for standard persistence. This adaptation is simplified by our observation that the incidence 
structure of a cell complex involved in level persistence computation can be derived easily from the given 
simplicial complex with a filtration. It is appropriate to mention that the level persistence considered here 
can be approached from the Zigzag persistence introduced in ||2l based on quiver representations. This and 
the Morse Novikov theory for circle valued maps from the perspective of quiver representations will be 
topics of further research. 

Note: After this work was done, we became aware that the "relevant numbers" for the level persistence 
of an M- valued map (as described later) can be derived from the invariants of the Zigzag persistence [21 . 
More recently we have noticed that they ai^e actually equivalent. Appendix D of a more recent paper Q 
makes this precise. 

'the general case can be reduced to almost integral case 

^values on integral 1-cycles are integer multiples of a fixed real number 
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2 l-cocycle 



Let X = {Xq, X) be a simplicial complex where Xq is the set of vertices, and is a set of finite subsets of 
Xq which satisfy the following properties : 

1. XqCX and 

2. Ifa e X then t C a ^ t e X. 

The subsets in X of cardinality A; + 1 are called /c-simplices. These subsets, considered as ordered k- 
tuples up to an even permutation, represent oriented A;-simplices. We denote sets of oriented A;-simplices by 
Xk- We continue to denote by X the canonical topological space (the geometric reahzation) associated with 
the combinatorial structure X. 

Broadly speaking, cochains are real valued functions that assign numbers to oriented simplices. We only 
consider 0- and 1-cochains which are sufficient for this exposition. A 0-cochain is a function restricted to 
vertices, / : A'o — > M. A 0-cochain can be identified with a continuous map / : X ^ M whose restriction to 
each simplex is linear (a linear map). The 0-cochain / is generic if f : Xq ^ Ris injective. 

The 1-cochains are are maps f : A*! ^ M whose domain is the set Xi of oriented edges of X. Recall that 
Xi can be regarded as the subset of Xq x Xq consisting of pairs {x, y) which ai^e vertices of a 1-simplex. 
The map f is a l-cocycle if it satisfies : 

1. f (x, y) = — f (y, x) for any ordered pair (x, y) € Xi, and 

2. if (x, y , z) e X2 then f (x, y) + f (y , + f (z, x) = 0; equivalently f (x, y) + ^{y,z) = ^{x,z). 
Clearly, a 0-cochain /, or equivalently a linear map, provides a l-cocycle 5f defined by 

8f{x,y) = f{y)-f{x). 

Any l-cocycle f represents a cohomology class < f >€ //^(X;R) and any such cohomology class is 
represented by a l-cocycle. Two 1-cocycles fi and £2 represent the same cohomology class iff fi — £2 = 5f 
for some 0-cochain /. 

An almost integral l-cocycle is a pair (f , a) where f is a l-cocycle whose values on integral 1-cycles 
are integer multiple of a fixed positive real a|j 

If St{x) denotes the star of the vertex x £ Xq (the star of any simplex is a sub complex), a l-cocycle f 
defines a unique function fx : St{x) — )■ R by the formulae fx{x) = and fx{y) = f{x, y) for any vertex 
y ^ xm St{x). Clearly, {fx — fy){z) is constant for any z in a connected component of St{x) n St{y). 

Thus, a l-cocycle can be thought as a collection of linear maps {fx ■ St{x) — t- M} for each vertex x, 
such that the difference fx — fy is constant on each connected component of St{x) n St{y). Therefore, one 
can regard f = {/a;, x G Aq} as a multivalued linear map. A l-cocycle £ is generic if all linear maps fx are 
generic, i.e., injective when restricted to vertices of St{x). 

2.1 1-cocycles and circle valued maps 

Consider a continuous circle valued map / : X ^ §1. Let p : E ^ = R/aZ be the map defined 
by p{t) = t(mod a), a a positive real number. For any simplex a £ X, the restriction /|cr admits liftings 
/ : (T — > M, i.e., / is a continuous map which satisfies p-/ = f\a- The map / : X — > is called linear if, for 

''in the language of algebraic topology the cohomology class < f >£ H^{X; R) has degree of rationality 1, or equivalently the 
image of the induced homomorphism Hi (X; Z) R is aZ C R. 
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any simplex a, at least one of the liftings (and then any other) is linear. Any linear map / : X ^ = M/aZ 
defines an almost integral 1-cocycle (f , a) by assigning to each edge e with vertices x, y, 

i{x,y) = fiy)-fix) 

where / is a lift of / restricted to e. One can show the opposite, that is, 

• an almost integral 1-cocycle {f : Xi ^ M, a) can be associated to a linear map / : X — )• = M/aZ. 

Covering associated with a cohomology class < f >G H^{X;R) : Regard X as a topological space. 
Choose a base point x G X and consider the space of continuous paths 7 : [0, 1] X with 7(0) = x, 
equipped with compact open topology. Make two continuous paths 71 and 72 equivalent iff the 71(1) = 
72(1) and the closed path 71 * satisfies < f > ([71 ★ 7^^]) = 0. Here denotes the concatenation 
of the paths 71 and defined by 72^^(t) = 72(1 — t), and [71 -k 72"^] denotes the homology class of 
71 -k 7^^. The quotient space X, whose underlying set is the set of equivalence classes of paths, is equipped 
with the canonical map vr : X — > X induced by assigning to 7 the point 7(1) G X. The map vr is a local 
homeomorphism and X is the total space of a principal covering with group G where 

G = img(< f >: Hi{X; Z)^R). 

When f is almost integral, G is isomorphic to Z. If X is equipped with a triangulation, then X gets a 
triangulation whose simplices, when viewed as subsets of X are homeomorphic by vr to simplices of X( 
when viewed as subsets of X). 

Construction of / and /. We construct the circle valued map / via its cyclic cover / as follows. 

Step 1. Consider vr : X — > X the principal Z— covering associated with the cohomology class < f > defined 
by f . This means that X is a simplicial complex equipped with a free simplicial action /i : Z x X — > X 
whose quotient space, X///, is the simplicial complex X. The map vr identifies to the map X — > X/^ = 
X and satisfies 7r*(< f >) = 0. 

Choose a vertex x of X and call it a base point. Notice that the vertices of X can be also described 
as equivalence classes of sequences {x = xq, xi, • • • xj\r-i,XN} with XjS being consecutive vertices 
of X (i.e. Xj, Xj+i are vertices of an edge). Two such sequences, {x = xq, xi, - ■ ■ , xat-i, x^v} and 
{x = yo,yi,--- yL-i,yL} ai-e equivalent if xn = yi and the sequence {x = zq, • • • , zn+l = x} 
with Zi = Xiifi < N and Zj+N = UL-j if 3 < satisfies 

^ i{zi,Zi+i) = 

0<i<L+N-l 

Step 2. Define the map / : — )■ M by f{y) := J2o<i<L~i ^{Vi^yi+i) where y G ^0 is the vertex coitc- 
sponding to the equivalent class of {x = ' ' ■ j Vl}- The description of X given above guarantees 
that / is well defined. Extend / to a linear map / : X ^ M. Observe that if yi and y2 satisfy 
T^iyi) = '^(^2) then f{yi) — f{y2) £ aZ. In addition if ei and 62 are two edges of X from yi to y[ 
and 2/2 to y'2 respectively with 7r(ei) = vr(e2), then f{y[) - /(y^) = f{yi) - f{y2). This imphes that 
if o"! and (T2 are two simplices with 7r(cri) = 7r(cJ2) = a and vr^s are the restrictions of vr to (Tj (vr are 
bijections on their image), then / • vrj~^ — / • vr^^ : u ^ M, is constant and this constant is an integer 
multiple of the fixed real number a. 
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Step 3. 



/ 



Observe that the map p ■ f : X ^ R (with p:R^S^ = R/aZ) 
factors through X//i = X inducing a map from X to S^. This 
is our linear map / whose associated 1-cocycle is f . The relation 
between / and / can be summarized by the commutative diagram 
on left. 



X 



p 



f 



3 Definitions of persistence 

This section defines persistence for almost integral 1-cocycles via circle valued maps. First, we expose the 
definitions for level persistence in analogy to standard persistence ||9l[l4|, and then extend them to circle 
valued maps. 

Let / : X — )• M be a continuous map. For simplicity suppose X is a nice compact space and / a nice 
map. This means that the homology Hr{-), r > 0, of all levels and sub levels are finitely generated. Define 



Persistence of a M-valued map. Standard persistent homology with coefficients in afield k |[T4l is the 
collection of vector spaces Br{t,t') := img(i/,.(X^t) — > ff,.(X|(/)), t < t' . It is understood that the 
homology Hr is considered with coefficients in the field k. The corresponding Betti numbers are /3r(i> t') '■= 
dimi?j.(i, t'). Equivalently, one can define persistent homology as the collection of vector spaces Kr{t) := 
Hr{X]_t) each of which is equipped with a filtration provided by subspaces: 



The corresponding Betti numbers ai^e Kr{t) := dimKr{t) and K,r{t,t') := dhjiKr{t,t'). The two type of 
informations are equivalent. Indeed, Br{t,t) = Kr{t) giving I3r{t,t) = Kr{t) and Br{t,t') is isomorphic to 
Kr{t)/Kr{t,t') implying = Kr{t) - Kr{t,t'). 

Definition 3.1 Let c e H.r(Xit) \ 0. One says that 

(i) c is bom at t', t' < t, if c e img {Hr{X^(^t'+e)) ^ Hr(X^t)) but c ^ img{Hr{Xj^(^t'-e)) ^ Hr(X^t)) 
for any positive e with t — t' > e > 0, 

(ii) cdiesatt", t" > t, if its image is nonzero in img {Hr{X^t) Hr{Xm/i_^'j)) but zero in img {Hr{X^t) 
Hr0^l{t" +t))) for any positive e with t" — t > e > 0. 

For any c G Hr{X^^t), let t~^{c) and t~{c) denote the death and birth times of c respectively where 
> t, t~{c) < t. For t' < t < t", one may introduce the numbers 

• fir{t', t, t") := the maximal number of linearly independent elements in HriX^t) so that t'^{x) = t" 
and t~{x) = t'. This also means that no linear combination of such elements is born before t', or dies 



The numbers firii', i,t") do not depend on t, hence can be written as firii', t")- The numbers /3r(i', i") 
determine and are determined by the numbers /ir-(t', t"), cf |l9l- The set of numbers f3r{t, i') or Krit), Kr{t, t') 
or t') are the relevant persistent numbers. 




r'((-oo,t]), 

f-\t), and 

rH[ti,t2]),t2>ti. 



Kr{t,t') := kei{HriX^t) ^ HriX^t')) with Krit,t') C Kr{t,t") C Kr{t) for t<t' < t" . 



before t" . 
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Simultaneous persistence. Suppose f± : X± [0, oo) are two continuous maps as above with (0) = 
/r^(0) = A. It will be useful to calculate the number 

• ujr{s, t) := the maximal number of linearly independent elements in Hj.{A) which die in X_ at s and 
in X_|_ at t (this means that no linear combination of such elements dies before s resp. t) is referred as 
the simultaneous persistence number of the pairs {/^ : — )• [0, oo)}. 

When X is a cell complex and / is a linear map restricted to each cell, there are effective algorithms to 
compute the numbers /ir(i', i")- These algorithms can be adapted to compute the number ujj.{s,t) for a pair 
: X^ — )• [0, oo) with X^ and A being cell complexes and being linear- on each cell as indicted in 
section 4. 

3.1 Level Persistence 

Level sets instead of sublevel sets define level persistence. Let 

L+(t;T) := ker(i7,(Xt) ^ F,(Xt,t+,)) and 
L-{t-T) := keiiHriXt) ^ HriXt^r,t))- 

Definition 3.2 The level persistent homology with coefficients in a field is the collection of vector spaces 
Lr{t) := Hr{Xt) equipped with two filtrations: 

L+{t;T) C L+{t;T') C L,.(t),r < t' and 
L~{t-T) C L~{t-T') C L,(t),T < t'. 

Consequently we have the relevant level persistence numbers: 

lr{t) :=dimLr(t), lf{t;T) := dimL^(i; r), and 
e^(t;r',r") := dim(L-(t; r) n r"))- 

Definition 3.3 Let c € H,.{Xt). One says that 

(i) c dies downward at t',t' < t, if its image is nonzero in img [HriXt) HriXtu^^t)) but is zero in 
img {Hr{Xt) Hr{Xt'^t)) for any e with t - t' > e > 0, 

(ii) c dies upward at t",t" > t, if its image is nonzero in img {Hr{Xt) — > Hr(Xt.t" -e)) but is zero in 
img {Hr{Xt) Hr{Xt^t")) for any e with t" - t > e > 0. 

For c € Hr{Xt), let t~^{c) and r~(c) be the upward and downward life time of c respectively, if c dies 
upward and downward at t + r"*" (c) and t — t~ (c) respectively. For r~ , r+ > 0, and t, we introduce the 
numbers : 

• r+) := the maximal number of linearly independent elements in Hr{Xt) so that t~^{c) = r"*", 

• T~) := the maximal number of linearly independent elements in Hr{Xt) so that t~(c) = r~. 

For tame maps, the numbers lf{t] r) determine and are determined by the numbers i'^{t; r). The numbers 
v^{t; t) can be calculated from standard persistence numbers ^r{', •) for another space Y equipped with a 
tame map : Y ^ M derived from (X, / : X ^ M). 

Similarly, the numbers er{t; t', t") can be derived from iOr{s, t) calculated for a pair {y^,g^) derived 
canonically from (X, /). Precisely, 
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1. t'jt (t; T~^) for map / is same as fj.r{t, t + r"*") for map g : Y — )• M where Y = oo and g = f\xt ^o- 

2. v~{t;T~) for map / is same as ij,r{—t,—t + r~) for map g : Y ^ R where Y = X_oo,t and 
9 = -/|x_oo,f 

3. er{t] T~, r+) is the same as the simultaneous persistence number uJr{T~ ,t~^) associated with Y+ = 

f-\[t, oo)), Y- = /-i((-oo, t]), g+ = /|y+ - t, and g^ = -f\y_ + t. 

3.2 Persistence for circle valued maps 

The standard (sub level) persistence can not be extended to circle valued maps because the notion of sub 
level is not well defined for circle valued maps. However, we can extend the level persistence to such maps. 
Let / : X — > = M/aZ be a continuous circle valued map. Consider 

• p : M ^ S^, defined by p{t) = t(mod a). 

• TT = vr-^ : X — )• X, the Z— principal covering associated to /, i.e., the pull back of the principal 
covering R ^ by the map /. Equivalently, it is the principal covering associated to the 1-cocycle 
defined by / as described in section 12.11 

• / : X — M is a map which satisfies p ■ f = / • vr-^. We make it unique by setting f{x) = for a 
chosen base point x of X. 

The map / is the infinite cyclic covering of / which would be same as the one constructed (section 12.1b 
from the 1-cocycle associated to /. Observe that 7r(Xt) = Xp(() = f^^{p{t)). 

Definition 3.4 Define L^{9; t) for f to be L^{t; r) for the map f where p{t) = 6. The definition involves 
the choice oft. It is easy to verify the independence of this choice. 

Now, we are ready to define persistence for circle valued maps. Let ^ G denote an angle and 
Xg = f~^{6). The level persistent homology with coefficients in a field for / is the collection of vector 
spaces Lr{9) := HriXe) equipped with two filtrations: 

• L+{e-T) C L+{e;T') C Lr{9),0 < r < r' < oo, 

• L-{e;T) C L-{e-T') C L^(0),O < r < r' < oo. 

Consequently we have the relevant persistence numbers: 

lr{e) := diuiLriO), 
1^{6;t) := dimL^(0;r), and 
e,(0;r',T") := dim(L-(0; r') n L+(0; r"). 

By the above definition, the numbers 1^{9;t) and er{0;T' ,t") for / are the numbers l^{t;T) and 
er{t; t', t") for /. Note that the numbers l^{6] r) and ej.{6; t't") can be computed from the level persis- 
tence numbers of / which in turn can be derived from persistence numbers for the maps g and simultaneous 
persistence numbers for the pair g^ as described in subsection 13.11 

3.3 Persistence for almost integral 1-cocycles 

Given an almost integral 1-cocycle (f, a), we associate the circle (S^ = M/aZ) valued map / : X — as 
described in section |2?T] and define the persistence of f as the level persistence of the circle valued map /. It 
is possible to shortcut the involvement of the map / and go directly to / to define persistence of f, but this 
might obscure the topological meaning of the definition. 
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3.4 Tame maps 

For finite calculations, one cannot allow the level sets change topology continuously. This is why we intro- 
duce the following restrictions of tameness on the maps. Similar conditions for standard persistence have 
been proposed before lH. In most practical situations, the tameness condition holds for the maps of interest. 

Definition 3.5 . A continuous map / : X — )■ Y, Y = M or §^ /i' called tame if there exists finitely many 
{ti, t2i • ■ ■ ^Af} that : 

(i) for any t ^ ti,t2, ■ ■ ■ t]\f there exists e > so that f : Xt-e,t+e — ^ [t — e,t + e] and the second factor 
projection Xt x [t — e,t + e] — )■ [t — e,t + e] are fiberwise homotopy equivalent, 

(ii) for any t (in particular for t = ti) there exists e > so that the canonical inclusion Xt C Xt-e,t+e 
a homotopy equivalence. 

The above definition can be considerably weakened by considering homology equivalence in place of 
homotopy equivalence, but in view of the fact that all our examples satisfy the above definition, we proceed 
with it. Generic smooth maps from a closed smooth manifold to R or are tame and so are piecewise 
linear maps from a simplicial complex to M or S^. If / : X — > is tame, then / : X — > M restricted to any 
Xt^ t2 is tame. 

Observe that, for a tame map / : X ^ M, the vector spaces Br{t,t') are completely determined by 
the vector spaces Br{i,j) := Br{ti,tj) and therefore (3r{t,t') and fir{t',t,t") by the numbers f3r{i,j) '■= 
I3r{ti,tj) and firihj, k) := fj,r{ti,tj,tk). The following relations between these numbers are well known H 

m. 

IJ.r{i, k) = firihj, k) = I3r{i, k — 1) — (3r{i, k) — I3r{i — 1, k — 1) + I3r{i — 1, k), 
^r{i,j) = ^ Atr(i',/) with dim Hr{Xit,) = Pr{hi)- 

j'>j, i'<i 

One may observe similar properties for level persistence. Let S2i = U and S2j-i be any number between 
ti-i and ti. Cleai^ly, the vector spaces Lr{t) and L^{t; r) are completely determined by the vector spaces: 
Lr{i) := Lr{si), L-^{i;k) = L+{si]Si+k - Si), L-{i;k) = L+{si]Si - Si_k)- Therefore, the numbers 
lf{t; r), i'^{t; r), and er(t; r", r"*") are determined by the numbers 
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— l^(ySi^SiJ^k Si) 




ii 


k) 


— Zj, (Sj, Si Si—k) 


yt 


{i 


k) 


= 1^+ [Si] Si+k - Si) 






k) 


= V-{Si\Si- Si_k) 


er{i 




k) 


— (sj, Si Si—j, Sj+fc 



Observe that 

l+{i;0) = l-{i;0) = 
/+(2i;l) = /~(2i;l) = 
er{i; 0, 0) = e^(i; 0, 1) = e^(i; 1,0) = er{2i; 1, 1) = 

and 

= X] ^rii'^k), l:;{i-j)= ^ v~{i;k). 

0<k<j 0<k<j 
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This shows that level persistence numbers can be computed from •) which in turn can be computed 
from fJ,r{',-) (section ISTI ). Similarly, the level persistence numbers er{i;j,k) can be calculated by the 
simultaneous persistence numbers ujr{si — Si-j, Si+k — Si) by the formula 

er{i;j,k)= ^ UJr{Si - Si_ji , Si+k' - Si) 

0<j'<j;0<k'<k 

4 Algorithm 

The algorithm to compute persistence for an almost integral 1-cocycle f follows the logical sequence that 
defines its relevant persistence numbers in the previous section. It considers the associated circle valued 
map /, and computes the level persistence for the cycMc covering /. We assume Z2-homology and describe 
how one can adapt a matrix version f6', '91 of the standard persistence algorithm |f8l| to compute level persis- 
tence for a M-valued map, then how to extend it to circle valued maps. In this development, we encounter 
subspaces of simplicial complexes which are cell complexes. We begin with its definition. 

Definition 4.1 A geometric cell complex X is a union of a collection X ofnon degenerate convex cells with 
disjoint interiors which satisfy the property that, for any cell a ^ X, all its faces belong to X. 

We consider finite cell complexes X equipped with a filtration J-^ := Xq C • • • X^ C X^+i • • • C X^ = X 
with Xj being sub complexes of X. 

Definition 4.2 A total order on X = {ai , . . . , an} is called topologically consistent if the condition A below 
is satisfied and filtration compatible if the condition B below is satisfied. 

• Condition A. cTj is a face ofaj implies i < j. 

• Condition B. ai G and CTj G X^./ \ X/^. with k < k' implies i < j. 

Given a filtration T, one can canonically modify any total order which satisfies Condition A into one which 
satisfies both conditions A and B. This is done in the following way: 

Find the first simplex a in the order which violates Condition B. Let r be the simplex immediately 
preceding a in the total order. The violation of condition B by cr implies that r G Xj\Xj_i and a G Xj\Xj_i 
where i < j. Permute a with r. Observe that condition A continues to hold; indeed the only possible 
violation of condition A is if r is a face of a which would imply that j < i. If a still violates Condition 
B in the new position, then continue moving a to left until it does not violate Condition B anymore. All 
other simplices which were initially preceding a do satisfy condition B. So, we have one less simplex which 
violates Condition B. 

Persistence algorithm lIH HI. The input to the algorithm is a matrix M that represents incidence struc- 
ture of a cell complex X equipped with a filtration. The algorithm derives the numbers fir{i,j) and the 
Betti numbers of X in case of Z2-homology. The numbers fir{i,j) determine the numbers f5r{i,j) for this 
filtration. Let X = {cri, . . . , Un} where ordering of the cells is topologically consistent and filtration com- 
patible. An entry M[i,j] in the incidence matrix M is set to the incidence number of crj and aj, that is, 
M[i,j] := I{ai, aj) where /(cjj, aj) is 1 if cjj is a face of codimension 1 of aj and otherwise. 

Let U be the class of 7i x 7i upper triangular matrices with zero on diagonal and all entries zero or 1 . The 
matrix M is in this class. Given a matrix M G U, for any column j, denote by low{j) the largest i so that 
M[i, j] 7^ 0. If the column j has all entries zero low{j) is not defined. 
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The persistence algorithm uses column additions to transform the matrix M within the class U with 
the additional property: no two columns j and / have low{j) = low{j') if they ai^e defined. A matrix 
with this property is called in "reduced form". The algorithm works by adding columns from left to right. 
Finally, the algorithm ends up with a matrix in the reduced form. The matrix in the reduced form obtained 
from M carries explicit homological information about X and the filtration. For example, the dimension 
of -ffr(X; Z2) is the cardinality of the set of zero columns j in the reduced form where aj has dimension r 
minus the number of nonzero columns / in the reduced form where aj' has dimension r + 1. To describe the 
number ^^ihi) we consider all pairs {ak,crk') with low{k') = k, and observe that /^^(z, j) is the number of 
such pairs which in addition satisfy 

1. k is an index corresponding to a cell in Xj \ Xj_i, 

2. k' is an index corresponding to a cell in Xj \ Xj_i, 

3. dimension of is exactly r. 



Now we describe how one may use the persistence algorithm to compute the numbers uJr{-,-) that we 
defined for simultaneous persistence. Suppose X is a cell complex and X='= C X are two subcomplexes with 
^ = X+ n X~. Suppose X+ and X' are equipped with finite filtrations CXf CX^ ■■■ Xt^ = X+ and 
X;: C X7 C X;7 • • • X" _ = X- where X^ = A. 



^2 "■N+ 

. . . y- _ V- „,Ua.-a v± _ A 
^0 — '^1 - '^2 "^N- 

One can put the cells of X in three groups I, II, III. The group I contains the cells of A, the group 
II the cells of X~ \ A and the group III the cells of X+ \ A. Suppose the cells are ordered according to 
condition A and condition B w.r.t. the filtration Xq C Xi C X2 • • • X7v-+Ar+ = X defined by Xj = X^^ for 
< z < A^" and Xj = X" U X^T^^^. for A^" + 1 < 7 < iV~ + A^+. The incidence matrix M is given by 

A B- C+ 

M = B with A being the incidence matrix for the complex A, the matrix M~ 
C 

being the incidence matrix for X~ , and M"*" 



A 


B- 





B 



, being the incidence matrix for X+ . 



are in reduced form. By running 
the relative reduced form can be 



k. 



A C7-+ 
C 

The matrix M is said to be in relative reduced form if M~ and M 
the persistence algorithm first on A and then continuing on M~ and M 
achieved. 

To compute the number we consider all triples {(7^,(7^1, a^") with low{k') = low{k") 

and observe that ujr{i,j) is the number of such triples which satisfy in addition 

1 . /c is an index corresponding to a r-cell in A 

2. k' is an index corresponding to a cell in Xj \ Xj_i, 

3. k" is an index corresponding to a cell in X^^^- \ Xj^jy--!- 

Next, we describe how we adapt the standard persistence algorithm to compute the relevant numbers for 
level persistence. First, we observe the following. If X is a simplicial complex, / : X — )■ M is a generic linear 
map, and to < ti, ■■■ ti <■■■ are the values of / on vertices of X, then the filtration of X by topological 
spaces /~^((— 00, tj]) is not a filtration by simplicial sub complexes. Apparently the above algorithm can 
not be applied. However it is possible to show that this topological filtration is homotopically the same 
(hence has the same persistence) as the filtration provided by the sub complexes Xj consisting of the union 
of simplices a so that /(a) C (— 00, tj]. 
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Computing level persistence for M-valued maps. In the case of level persistence, the calculation of 
the numbers and u^{i,j) for a tame map is reduced to the calculation of the numbers •) for 

subspaces of X, precisely X^, Xj oo, X_oo,tequipped with the appropriate function , and X^j indicated at 
the end of section [3?T] In what follows we describe how to derive the matrix M for X^, X^.oo ^-oo.t, "^tify 
from the matrix M for X. 

Let X be a simplicial complex and / : X — )• R be a generic linear map. Observe that X^, X^ oo ^-oo,t, 
X^^ ^2 cell complexes as indicated in the previous section and not simplicial complexes. However, the 
incidence structure of their cells can be described in terms of the incidences of simplices in X. This key 
observation allows us to construct the incidence matrix of the cell complexes X^, X^ oo "^-00,11 ^ti,t2 from 
the matrix M{X) as detailed below. For any simplex a, let [a] denote the closed interval which is the convex 
hull of the numbers {f{v), v a vertex of a}. 

Cell complex Xf For any t introduce the matrix M{t) as the minor of M(X) consisting of the rows and 
columns i with t € int [cji] and regard the simplex as the cell ai with dim (ctj) = dim (cjj) — 1. If there 
are ^ < n such simplices, M{t) is an £ x £ matrix. Clearly, M{t) respects the order of the cells induced by 
the order of simplices used for M (X). 

• If no vertex takes the value t, then M{Xt) = M{t) is the incidence matrix of the cell complex Xf . 

• If there exists a vertex v (unique since / is generic) so that f{v) = t, then the incidence matrix M (Xt) 
is an (£+ 1) X (£ + 1) matrix with one additional row and column corresponding to the vertex v viewed 
as a cell of dimension which should be indexed before all other cells. The entry for the pair {v, aj) 
in M (Xt) = 1 if cjj is a 2— simplex which has t; as a vertex and otherwise. 



Cell complexes Xt^oo, ^-oo,t, ^ti,t2- To describe the incidence matrices which correspond to the cell com- 
plexes Xf oo> X_oo,t, and Xf^ fj with ti < t2, consider the collection of simplices: 



X{-(x>,t) : = 



= {a, [a] n (t, 00) ^ 0} 
{a, Hn(-oo,t)/0} 



Let the matrices M(t, 00), M{—oo, t), and M{ti , be the minors of M (X) whose rows /columns are 
the simphces in the sets X{t, 00), X{—oo, t), and X{ti,t2) respectively. We keep the dimension of the cell 
same as that in the matrix M(X). The matrices M(Xf 00) and M(X_oo,t) are matrices of the form 



M{t) 




M{t, 00] 



and 



M{t) 




M{-oo,t) 



respectively , 



The matrix B'^ is a s x m matrix where Xf has s cells and X{t, 00) has m simplices. We set B^[i,j] = 1 
if the i-th cell in Xf and the j-th simplex in X{t, 00) ai^e either a and a respectively, or a vertex of X and a 
1-simplex incident to it respectively. All other entries in i?"^ are 0. Define B^ analogously by considering 
simplices in X{—oo, t). The incidences among the cells in Xf 00 are of three types, the ones among the cells 
in Xf, the ones among the cells in Xf 00 that are not in Xf, and the ones among the cells one of which is in 
Xf and the other is in Xf 00 but not in Xf . The matrices M{t), B^ and M{t, 00) capture these three types of 
incidences respectively. Clearly, similar observations can be made about the incidences among the cells in 



Let Af(ti,t2) denote the minor of MiX) whose rows and columns correspond to the simplices in 
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X{ti,t2). The matrix M(Xt, ,42) is 
















B~ 

M{ti,t2) 



where the matrix has the entries as defined in the previous case with X^^ and X{ti,t2) playing the roles 
of Xt and X{t, 00) respectively. The matrix B~ is analogous. 

The order of the cells as suggested by the matrices M{Xt), M{X-oa,t), M{Xt,oo) and M{Xti,t2) satisfies 
Condition A but not necessarily Condition B w.r.t. the filtration induced from the map g (end of section IXTI ). 
If necessary it can be canonically changed as indicated at the beginning of this section. 

Computing level persistence for a circled value map. Let / : X — be a generic linear map. For 
simplicity in writing we suppose that f{a) is an angular interval smaller than vr. We need the incidence 
structure of the simplices in the covering space X and the linear R-valued map /. However, the space X is 
infinite. So, we compute only a finite subspace of X which is sufficient for computing relevant persistence 
numbers and the restriction of / to these subspaces. 

We show how to construct a simplicial complex from X that contains Xt,t4-27rfc- The constructions of 
Xt_27rfc,t, Xt_27rfci,t+27rfe2 are analogous. 

Let X denote the set of simplices in X. For any G decompose Af as a disjoint union X = 

UL^ U d^L^ U d^L^ where (see Figure □ 

• CP consists of the set of all simplices whose closure do intersect the level Xq . Let be the simplicial 
complex generated by simplices in CP ^ 

• is the set of simplices which do not belong to L^. Let denote the simplicial complex generated 
by the the simplices in and consider T^ n L^. This simplicial complex is the disjoint union of two 
simplicial complexes d^VP and characterized by /(cj) < for a G d-VP and /(cj) > Q for 



• d±CP represent the simplices in d±\P . 

Our purpose is to build a collection of simplicial complexes which are equivalent with the space Xt^t+2TTk 
where p{t) = 0. 



Figure 1: Complex X with level Xq on left. Complexes X (1), X (2), • • • providing a filtration of the total 
space X on right. 

Introduce a nested sequence of simplices ^^(0) ^ ^^(1) C • • • X^{k) as follows. Since we will repeat 
copies of each of the sets , d.C^, and d+£^, let T\n), C\n), d_C%n), and d+C^n) denote their 

n-th copies respectively. Taking C^{0) = C\ T^(0) = d^C^Q) = d+C^, define inductively, 

X\0) = d^C^ 
X^{n+1) = X\n)UC\n)Ud+C^{n)uT%n)Ud^C^{n + l) 
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Taking Io{a, r) = I {a, r) and assuming r to be a face of a of codimension 1, the incidences among the 
simpUces are described by 



/„+l(cr,r) 
/„+i(fj,r) 
/„+i(o-,t) 



I„((j,r) ifae X^{n) 

I {a, r) if cj € C^{n) U 9+£^(n) viewed as simplices of 

/(cr, r) if cr € T^{n) U d^C^{n + 1) viewed as simplices of T^. 



In all other cases /„+i(cj, r) = 0. Notice that each forms a simplicial complex X^(i) (Figure [T]). To 

describe / it suffices to provide its values on vertices. We write 



and ■Po for the subset of vertices in V. We write VQ{n) for the n— th copy of Vn and define /(n) : Vo{n) — )• M 
by f{n) := / + 27rn where / = • / with p : [t,t + 27r) — whichQ sends t to 9. Once defined on 
vertices, / is extended by linearity to each simplex of the simplicial complex X^(n). 

Note that an order of the simplices of X satisfying condition A induces an order on the simplices of 
T^(n), d±C^{n) and C^{n) and by juxtaposition an order on the simplices of X^(n) which continue to sat- 
isfy condition A. It implies that one can build a matrix M(X^(?i)) which satisfies condition A by juxtaposing 
the minors of M(X) that represent , d±L^, , and their copies in an appropriate order. Note also that 
X^(77,) is a sub complex of X and therefore the restriction of / provides tame maps on each of these spaces. 
The columns and rows of M {X^{n)) can be reordered so that they become filtration compatible with f or g 
by the method indicated at the beginning of this section. 

Algorithm for almost integral 1-cocycIes. If (f , a) is an almost integral 1-cocycle, hence given by M(X) 
and a real number f (e) for each edge, we have to provide the map / : X — hence a matrix M (X) and 
an angle valued map on vertices and calculate the level persistence. For this purpose we choose a base point 
vertex x and assign to it the angle value 0. For any other vertex y choose a sequence of consecutive vertices 
{x = yo,yi, - ■ ■ yL~i,yL = y} and assign to y the angle (Eo<i<L-i ^ yi+i)) ^od a)2Tr/a. The result 
is independent of the choice of x. We continue then as described above. The total number of copies of L 
and T involved in the calculation of aU persistence numbers described here is less than the cardinality of X. 



"^p is bijective and continuous, but p ^ is not continuous 



V = d^r uru d+u u V 



■e 
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Appendix 



In this Appendix we show that, for a tame map / : X ^> M, the relevant level persistence numbers determine 
the persistence numbers for both / and — /. We suppose X is compact so / is bounded from above and 
below. 

Let {ti,t2, ■ ■ ■ tk} be the critical values of / and let {t'l, ■ ■ ■ i^+i} regular values so that t[ < ti < 
t'2 < t2 < ■ ■ ■ t'l. < ■ ■ ■ tk < ifc+i- Consider S2i = ti and S2i-i = t'-. Clearly = and Xsj^+i = X, 
therefore Xsi,^. = X_oo ,;. and X^. S2fe+i = X^.^oo- Note that the tameness of / implies that the inclusions 
X^. C 'Xs^,Si+i for i even and X^^^^ C 'Xsi,s,+i for i odd induce isomorphism in homology. 

The relevant level persistence numbers are lr{si), lf{si] ■iz{si±k — Si) and e,.(sj; Si — Si-k', Si^k" — Si) 
while the relevant persistence numbers for / are Kr{i) = dim Hr{'X-oo,Si) and 
K,r{si, Si+k) = dimker(ii'r(X_oo,sJ i^r(X_oo,s,+ J) while for -/ are dimHr{X^s„oo) and 
dimker(iJr(X_s,,oo) HriX_s,+k,oo))- 

To show that the numbers lr,lf,er determine the numbers we will use the algebraic proposition 
below: 

Proposition 4.3 Let 

^ yl^ ^ i?^ ^ . . . ^„ ^ i?„ ^ C„ ^ An-i . . ^ ^0 ^ i?o ^ Co (*) 

< n < N be an exact sequence of vector spaces. Then any three independent collections of num- 
bers {dim An}, {dimS„}, {dimC„}, {dim(ker a^)}, {dim(ker /3„)}, {dim(ker 5^)} determine the other 
three. (Here independent means not constrained by the obvious equalities (1), (2), (3) below. 

We have the following two long exact sequences: 

1. The Mayer- Vietoris sequence associated to Xt^t" = Xj^^/ U Xf^t", ^t,t' H Xti^t" =Xt',t<t'< t". 

■ ■ ■ Hn{Xt^t')®Hn{Xti ^tii) ^ HniXt,t") Hn-l{Xti) Hn-l{Xt^t')®Hn-l{Xti ^tii) ■ ■ ■ . 

2. The long exact sequences of the pairs Xt j/ C Xt^t" and Xf^t" C Xt^t", t <t' < i" ■ 

{I) • • • Hn{Xt^t') Hn{Xt^t") ^ Hn{Xt^t//,Xt^t/) ^ Hn-l{Xt^t') ' ' ' ■ 

{II) ■ ■ ■ — Iln0^t',t") Iln0^t,t") ^ Iln0^t,t" ,^',1") Hn-lO^t' ,t") ' ' ' • 

For t <t' < t", we also have by excision theorem 

(///) Hn iXt,t" ,^t,t') = Hn (Xt' , Xj/ ) , Hn iXt,t" , ^t' ,t") = IIn{'^t,t' , Xf ) • 

Proof: [Proposition 14.31 : Note first that the exactness (*) implies the following equalities 

(1) dim An = dim ker q„ + dim ker I3n 

(2) dim Bn = dim ker /3„ + dim ker (5„ (1) 

(3) dim Cn = dim ker 6n + dim ker On-i 

Then we obtain 
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• dimj4„,dimi?„,dimker a^, < n < in view of (1), (2) ,(3) determine dim ker /?„ , dim ker 5^ , 
and dim Cn <n < N. 

• dim^n, dim dimker /3„, < n < in view of (2), (1), (3) determine dimker(5„, dimkera„ 
and dim C„ ,0 < n < A^. 

• dim^n, dim C„, dimker a„, < n < in view of (1), (3), (2) determine dimker /?„, dimker(5„, 
and dim Bn ,0 <n < N. 

• dim^n, dimC„, dimker /3„, < n < A^ in view of (1), (3), (2) determine , dimker a„, dimker 5„ 
and dim Bn ,0 <n < N. 

• dimker Q„, dimker /3„, dimi?„, < n < A^ in view of (1), (2), (3)determine dim A^, dimker (5„ 
and dim C„ ,0 < n < A^. 

• dimker Q„, dimker dim i3„,, < n < A^ in view of (3), (3), (2) determine dim C„,, dimker /3„ 
and dim An ,0 < n < N. 

• dimker a„, dimker dim < n < A^ in view of (1), (2) , (3)determine ker dim/5„, dimi3„ 
and dim C^^ ,0 < n < A^. 

• dimker a„, dimker /3„ dimker 5„, < n < A^ in view of (1), (2) , (3) determine dimA„, dimi?„, 
and dimC, < n < A^ 

• dim An, dim ker dim ker 5„ ,0 < n < A^ in view of (1), (2) , (3) determine dim ker /3„, dim Bn 
and dimC„, < n < A^. 

• dim A„, dimi?„, dimC^, < n < A^ in view of (1), (2) , (3) and the fact that dimkeroAr = 
determine dim ker a^, dimkev Pn, dim ker (5„, < n < N. 

All other possible situations can be recovered from these cases. ■ 
To conclude that the level persistence numbers determine the persistence numbers we proceed as follows. 
Use the equality dimffr(Xs.^si+i) = dim//r(XsJ if i even and dim Hr{Xsi,Si+i) = dim. Hr{Xs^_^_j) if i odd 
which follow from tameness, and apply Proposition 14.31 to the exact sequence (1). In this way we derive 
from the numbers lr{si) and er{si; Si — Si-w, Si+k" — Si) the numbers dimHr{Xs^^s^+^.) for all i, k. 

Apply Proposition 14. 3| to the long exact sequence 2(1) for t = t', (resp. 2(11) for t' = t"), and in view of 
(III) derive from the numbers lr{si) and l^{si; Sj+fc — Sj) (resp. l~{si] si — Si^k)) the numbers 
dimHr{Xs,,s,+,+k,y^s,+,.,s,+,+u) (resp. dimF,.(X,^,,^^^^,, X,^,,.+J). 

Apply Proposition 14.31 to the long exact sequence 2(1) (resp. 2(11)) and derive from dim Hj-i^s^^Sij^,,) 

dimHr{X,^^Si+r+k,y^s,,s^+^ (resp. dimi7^(X,^,,^_^^_^,,X,^^^,,._^^^J) the numbers dimker(F^(X,^,,,^J ^ 
^r(Xs„s,+,+J) (resp. dimker(i7^(X,.^^,,.^^^J ^ F^(X,^,,^^^^J)). Taking Si = si (resp. Si+k+r = 
S2k+i) one obtains the persistence numbers for / and — /. 
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